A numerical method for detecting incommensurate correlations in the Heisenberg 

zigzag ladder 
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We study two Heisenberg spin-1/2 chains coupled by a 
frustrating "zigzag" interaction. We are particularly inter- 
ested in the regime of weak interchain coupling, which is 
difficult to analyse by either numerical or analytical meth- 
ods. Previous density matrix renormalisation group (DMRG) 
studies of the isotropic model with open boundary conditions 
and sizeable interchain coupling have established the presence 
of incommensurate correlations and of a spectral gap. By us- 
ing twisted boundary conditions with arbitrary twist angle, 
we are able to determine the incommensurabilities both in the 
isotropic case and in the presence of an exchange anisotropy 
by means of exact diagonalisation of relatively short finite 
chains of up to 24 sites. Using twisted boundary conditions 
results in a very smooth dependence of the incommensurabil- 
ities on system size, which makes the extrapolation to infinite 
systems significantly easier than for open or periodic chains. 



I. INTRODUCTION 

In recent years several frustrated quasi one dimensional 
magnetic compounds have been identified and studied 
experimentally In the one-dimensional phase, i.e. 

for temperatures above the magnetic ordering transition, 
frustration is expected to lead to incommensurate cor- 
relations. Precisely how frustration gives rise to incom- 
mensurabilitites for the extreme "quantum" case of spin 
S = 1/2 is at present not well understood. 

A paradigm of a frustrated quantum magnet is the 
spin-1/2 Heisenberg antifcrromagnetic chain with near- 
est neighbour exchange J\ and next-nearest neighbour 
exchange J2. This model is equivalent to a two- leg lad- 
der (see Fig. 1), where the coupling along (between) the 
legs of the ladder is equal to J 2 (Ji). 

The Hamiltonian is given by: 

H = £>i(SfS? +1 + SfS? +1 AS?S? +1 ) 

i 

+ ]T J 2 (Sf Sf +2 + SfSf +2 + AS*S? +2 )\, (1) 

i 

where we have allowed for an exchange anisotropy A. 
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FIG. 1. Schematic representation of a zig-zag ladder (top) 
and its equivalent chain with first and second nearest-neighbor 
exchange interactions (bottom). 

The zigzag ladder model is believed to describe the 
quantum magnet SrCu02 above the magnetic or- 
dering transition. The ratio of exchange constants is 
estimated as \J\/Ji\ ~ 0.1 — 0.2 so that the inter- 
chain coupling is significantly smaller than the exchange 
along the legs of the ladder. A second well-studied ma- 
terial with zigzag structure is CS2CUCI4 [0]. However, in 
CS2CUCI4 the chains appear to be coupled in an entire 
plane and no pronounced ladder structure is present. 

The model ([!]) in the regime J2 > J\ has been studied 
previously using both numerical [0J^j and field-theoretical 
Ig-Q methods. The main DMRG results for the spin 
rotationally symmetric case (A = 1) and not too large 
values of J%jJ\ are the following M 

• The ground state is doubly degenerate and is char- 
acterized by a nonzero dimerisation d — (S^n ■ 

• The equal-time correlation function (S n ■ S±) ex- 
hibits an oscillating exponential decay at large spa- 
tial separations. The characteristic angle associ- 
ated with these oscillations, i.e. the incommensu- 
rability, is related to the correlation length £ by 
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This connection of incommensurability and corre- 
lation length also fits into the picture emerging 
from the renormalisation group analysis (valid in 
the limit J2 3> Ji) of ||, which yields the simul- 
taneous divergence, with fixed ratio, of the related 
coupling constants. 

The regime J2 3> J\ is very difficult to analyse numer- 
ically as both the dimerisation and the incommensurabil- 
ity become very small. In particular, the DMRG analysis 
of H did not consider this regime. At the same time the 
field theory studies of || suggest that for J 2 ^> J\ and 
A I < la different type of physics may emerge: there are 
still incommensurate correlations and dimerisation, but 
in addition there is also "chiral" order 

(^2n^2n+2 — ^2n^2n+2) 7^ I 

{S2 n S 2n -l — ^2n^2n-l) 7^ ■ (3) 

Such type of order is forbidden in the isotropic case 
A = 1. It clearly would be interesting to numerically 
analyse whether or not such a new phase indeed exists. 
Given the difficulties, mainly due to finite-size effects, 
in accessing the relevant parameter regime by numeri- 
cal methods, we propose in essence to utilise finite-size 
effects to extract information on the incommensurabil- 
ity present in the system. This is done by studying the 
effects of twisted boundary conditions on the energy lev- 
els. Our main purpose is to establish the viability of 
our method by carrying out exact diagonalisation of fi- 
nite clusters of up to 24 sites. In order to obtain definitive 
results on the zigzag chain in the most interesting param- 
eter regime, larger systems need to be studied, possibly 
by implementing TBA into a DMRG algorithm. 

II. THE METHOD 

We found the ground state of the system by Lanczos 
diagonalization of rings of size L = 12, 16, 20 and 24 sites 
in the subspaces of total spin projection S z — 0, 1, for 
each total wave number K. We used twisted boundary 
conditions (TBC) S~[ +L = e^Sf. This means that each 
time a spin down traverses one particular link, it acquires 
a phase e l * (e~ 1 *) if it moves to the right (left). In the 
fcrmionic representation, after a Wigner-Jordan transfor- 
mation = c\ exp(nr X)j<i c ] c j)i and a gauge transfor- 
mation c\ — e l ®/ L fl , the problem becomes equivalent to 
a system of fermions (/) on a ring threaded by a flux 
where is the number of spins down. The advan- 
tage of the TBC for our purposes is that the allowed total 
wave vectors are: 

Kn ^) = ^n+jN l , (4) 

with n integer. Thus, varying we have access to a con- 
tinuum of possible wave vectors, even if we are working 



on a finite system f|. As an example, the ground state 
of Eq. ([j]) for N± = 1 is known exactly. In the thermo- 
dynamic limit, for 1/4 < a = J%l J\ < 1/2, the ground 
state is two fold degenerate with incommensurate wave 
vectors K = ± arccos[— l/(4a)] (6). For small systems, 
these wave vectors are only accessible for discrete partic- 
ular values of a if PBC are used. Instead, for any a, the 
exact JsT's and ground state energies are reproduced, if 
the energy of a small ring is minimized as a function of 
flux <f> and discrete wave number (n). 

For |A| ^ 1, the transverse spin correlations dominate 
at large distances. We denote the ground state of the 
system by \g). It lies in the S z — sector and its wave 
vector for PBC (or TBC if the energy is minimized over 
$) is always Kq = or tt. For the transverse spin corre- 
lations we can write: 

G{1) = (g\S+Sr +l \ g ) = I £ e-* l \(e\S-\g)\*, (5) 

e,q 

where S~ is the Fourier transform of Sf, and the sum 
over |e) runs over all excited states. 

By symmetry, for each q, only excited states with 
S z = — 1 and K\ = Kq + q have nonvanishing matrix 
elements in Eq. (||). We now assume that the large- 
distance asymptotics of C(l) is determined by the lowest 
excited state in the S z = — 1 sector. The difference in 
wave number between this state and the ground state 
then gives the incommensurability g m i n . 

For a finite system, in order to represent a continuum 
of wave vectors, we minimize the energy in the S z — 
— 1 sector E\ (Ki , <I>) , as a function of flux and allowed 
discrete wave vectors for this flux (see Eq. (ji))). The 
wave vector of the excitation is taken as: 

9mi„ = = #l($min) ~ #o($min), (6) 

where $ m i n is the flux which minimizes Ei(Ki,<f>), and 
Ki(&) is the wave vector of the state of lowest energy in 
the S z = — 1 sector at flux (f>. The wave vector g m i n gives 
the period of the oscillations of the spin correlations and 
corresponds to the pitch angle 9 of the classical spiral 
density wave Q. 

We have also investigated the energy gap of the spec- 
trum. For a finite system we evaluate it as: 

A 9 = E^Ki, $ min ) - E Q (K , $; nin ). (7) 

Here ^' min is the flux which minimizes the ground-state 
energy Eq. There are alternative expressions to q m in and 
A g which converge to the same value in the thermody- 
namic limit. Our experience suggests that the ones we 
chose have the fastest convergence. 

To test the method, we have studied a dimerized half 
filled system of non interacting spinless fermions. This 
model is the fermionic version of an XY model plus ad- 
ditional interactions: 
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= + H.C.) 

i=l i 

+ ^E(- 1 ) J ( c li^+ H - c -)- (8) 

i 

We have taken the parameters t\ = 1, ti = 0.4, £3 = 0, 
<4 = —0.2, V = 0.3, in such a way that the upper (empty) 
band has minima at incommensurate wave vectors q min = 
±0.2356 7r, and the lower (full) band has its maximum 
at q = 0. There is an indirect gap A g = 0.12477. 
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FIG. 2. Incommensurate wave vector (top) and gap (bot- 
tom) as a function of system size, obtained using periodic and 
twisted boundary conditions, for the toy model Eq. (8). 

In Fig. 2, we compare the results for ^min and A g 
in finite systems of L sites, with L/4 integer, and 12 ^ 
L ^ 64, between PBC and TBC. The gap is calculated as 
A g = Li($i) + L_i($_i)-2L ($o), where E^) is the 
ground-state energy for i added particles. For TBC the 
fluxes <&i are those which minimize while for PBC 

= 0. Since q min is near 7f/4, and the latter is one of the 
allowed wave vectors of PBC for L multiple of 8, small 
periodic systems with L/8 integer have the minimum en- 
ergy for one added particle, when this particle has wave 
vector ±7r/4. As long as 2ir/L is larger than 7r/4 — \q m i n \ 
(small systems), the results for PBC are better if L is 
multiple of 8. However, the oscillations with increasing 
L for PBC make a finite-size scaling difficult. Although 



q m i n also oscillates with L for TBC, the oscillations are 
smaller and the convergence to the thermodynamic limit 
is much faster. Using TBC, for L ^ 52, the error in q m i n 
and A g are below 0.01%. For the maximum size of the 
system used in our Lanczos diagonalization of Eq. (Q) 
(L = 24) , the error in q mul is ~ 2% and that of A g is of 
the order of 10%. 



III. RESULTS 

1. a) Isotropic case 

The incommensurate spin correlations and spin gap for 
A = 1 and J%j J\ < 3 have been studied previously by 
DMRG @,§. In this section, we compare our results with 
these ones, and extend the study to 3 < J2/J1 < 30, a 
region which is very difficult to reach with other methods. 

We have used a linear extrapolation in 1/L of the data 
for the angle 6 = q min for L = 12, 16, 20, and 24. We 
have chosen L/4 to be an integer in order to avoid frus- 
tration of antiferromagnetic interactions for large J2. A 
quadratic fit gives smaller values of 9, which underesti- 
mate the DMRG results. The comparison of available 
DMRG results and those obtained using TBC as de- 
scribed in the previous section is included in Fig. 3. For 
J2/J1 < 0.7 we do not obtain any incommensurability. 
This is probably a finite size effect, since for J2/J1 = 0.7, 
we obtain 9 — ir for L = 12, but incommensurate val- 
ues of 9 for L > 12. We have disregarded the value for 
L = 12 in the extrapolation when J2/L1 = 0.7. 
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FIG. 3. Incommensurate angle as a function of J2/J1 for 
A = 1, obtained using twisted boundary conditions. The 
DMRG results of Refs. 7,8 are also shown. 

For J2/J1 < 1 and J2/J1 > 1.8 our results are in bet- 
ter agreement with those of White and Affleck || than 
with those of Bursill et al JtJ. In the remaining region 
both DMRG results are very similar. In general, the dif- 
ference between the nearest of the results of the three 
calculations for 9 — 90° is of the order of 20%. For exam- 
ple for J2/J1 — 0.7 our results and those of Ref. are 
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respectively 31°, 28° and 21°. For J 2 /Ji = 2 the corre- 
sponding values are 1.00°, 1.19°, and 1.36° . In general, 
comparison with DMRG results and the difference be- 
tween different extrapolation methods suggest that the 
error in 9 — 90° using TBC is roughly of the order of 
20%. In view of the simplicity of our method compared 
to DMRG calculations, we believe that our results are 
satisfactory. Note that to detect an incommensurability 
of 1° without using TBC, the size of the system should 
be of the order of L ~ 360 ! 



We have also used (0) to study the gap A g . Due to 
the smallness of the gap, the finite-size effects for the 
relatively small system sizes we consider are too large to 
allow us to obtain reliable values for A g . Fig. 5 shows 
the size dependence of A g . From the difference between 
linear and quadratic extrapolation, we estimate the error 
in the extrapolated gap to be of the order of O.lJi, while 
for any value of J2/J1, A ff < 0.5 J\, Within this error, our 
results agree with those reported by White and Affleck 
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2. b) Anisotropic case 

In H the zigzag ladder was studied by means of a field 
theory approach in the regime J2 S> J\ ■ A mechanism for 
generating incommensurabilities was identified and anal- 
ysed quantitatively for the case of two coupled XX chains 
(A = 0). It was found that spin correlations exhibit a 
very slow power-law decay and are incommensurate 
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FIG. 4. Same as Fig. 2 for J2/J1 > 2 in logarithmic scale. 
Dashed line is the function e - 2 -0 05J 2 / Ji 

In Fig. 4, we show the angle as a function of J2/J1, 
for 2 < J2/J1 < 30 on a logarithmic scale. In this region 
the incommensurability is very small, and therefore, as 
explained above, very hard to obtain by alternative nu- 
merical methods. For large J2/J1, the deviation of the 
angle from n/2 is expected to be of the form MM 



-90° = cxp(-a- 6J2M) 



(9) 



A linear fit of ln(9 — 90°) as a function of J2/J1 in the 
interval [15,30] gives within a few percent a = 2, and 
6=1/20. 
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FIG. 5. Size dependence of the gap for A = 1 and several 
values of J%l 3\. Dashed (full) lines are quadratic (linear) fits 
of the data. 



where j — 1,2 and the deviation of the pitch angle from 
ir is k oc {Ji/h) 2 - The analysis of || also implies the 
existence of local magnetisation currents around the ele- 
mentary triangular plaquettes of the ladder. The findings 
of H were questioned in Jl3|], where the squares of the 
local magnetisation currents were computed numerically 
and found to decrease with system size for open chains 
of up to 16 sites. 




FIG. 6. Incommensurate angle as a function of A for 
J2/J1 = 10. 

In an attempt to resolve this controversy we have set 
J2/J1 = 10 and studied the variation of the incommen- 
surate angle with the anisotropy parameter A. As shown 
in Fig. 6, we find that 6 — 90° increases considerably as 
A is decreased from the isotropic case A = 1. This is in 
agreement with the field-theory prediction of ||. 

We futhermore have determined the dependence of the 
incommensurability on J2 / J\ ■ The results are shown in 
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Fig. 7. For large values of J2/J1 our numerical results 
are well fitted by 



90° = 2.785° ^ 
J2 



(11) 



This is in disagreement with the prediction of |j) . 

The disagreement between the predictions of H and 
the numerical results ( |lT| ) and pl| could either be due to 
a defect in the mean-field solution of Q or be an artifact 
of the limited system sizes used in the numerical com- 
putations. In fact, the analysis of Q predicts a gapless 
phase at A = so that it is conceivable that numerical 
results for small clusters are plagued by finite-size effects. 



CD 

CD 0.25 

O) 
CD 

T3 0.20 
o 

S 015 
CD 

0.10 
0.05 h 
0.00 



10 15 20 25 30 

J 2 /J, 

FIG. 7. Incommensurate angle as a function of J2/J1 for 
A = 0. Dashed line is the function 2.785° J1/J2. 

There is indeed evidence suggesting that finite-size ef- 
fects are still significant for L = 24. We find that the 
ground state in the S z = 1 sector has a lower energy 
than the first excited state in the S z = sector for L = 24 
and J2~2Ji. On the other hand, the DMRG studies of 
H show that in the isotropic case and for long lattices 
the two lowest levels are both in the S z = sector (de- 
generate ground states corresponding to different signs of 
the dimerisation) . 

We note that the presence of such finite-size effects 
does not necessary imply that the extrapolated results 
for the incommensurability are incorrect. In order to 
resolve this issue it is necessary to study significantly 
longer lattices. 



IV. SUMMARY AND DISCUSSION 

We have calculated the incommensurate wave num- 
ber in the next-nearest-neighbor Heisenberg model with 
anisotropy A, by exact diagonalization of rings of up to 
24 sites. We have used twisted boundary conditions and 
assumed that the incommensurate spin fluctuations for 
A I ^ 1 are determined by the lowest excited state for 
total spin projection S z = ±1. 

The method is able to detect incommensurate angles 
9 ~ 0.03°. This corresponds to a wave length of the order 



of 10 000 sites and is impossible to detect by using alter- 
native numerical methods. However, for certain param- 
eters (J2/J1 < 0.7 in the model), our method is unable 
to detect the incommensurability although it is rather 
large. On the other hand, the method does not predict 
incommensurabilities in cases where it is known that none 
exist. Also, in general, the extrapolated vales of 9 seem 
to be underestimations, as compared with known DMRG 
results. This is also the case for the toy model (Eq. (||)), 
represented in Fig. 2, where a linear extrapolation gives 
an underestimation of <7 m i n by ~ 30% if the results are 
limited to 24 sites. 

The advantage of using TBC for facilitating a finite- 
size scaling analysis has been noted previously in e.g. Q , 
but their use for detecting incommensurabilities is to the 
best of our knowledge novel. 

In spite of the limitations of the size of the cluster, the 
values of 9 obtained with our method are in reasonable 
agreement with known DMRG results in the isotropic 
case. For this case, we have also studied the region 
J2I J\ > 3, which is very difficult to reach by alternative 
methods. For sufficiently large J2/J1, 9 — tt/2 decays as 
cxp(— b J2/J1) as predicted by field theory We obtain 
that the constant b ~ 1/20. 

In the anisotropic case, we obtain that 9 increases with 
decreasing A, in agreement with Ref. ||. However, we 
find a linear dependence of 9 — n/2 on Ji / J 2 for J 2 ^> J\ , 
in contrast to the quadratic behaviour predicted in ||. 
This discrepancy may be due to finite-size effects. 

In summary, we have shown that incommensurabilities 
can be detected by diagonalizing finite-size clusters with 
TBC. The main advantage of the method is that the 
dependence of the incommensurability on system size is 
very smooth and allows extrapolation from results for 
relatively short chains. 

It would be very interesting to implement our TBC 
method in a DMRG algorithm and study the anisotropic 
zigzag chain for much larger sizes. 
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